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The in-plane resistivity of the high-temperature oxide superconductor La2-c C Sr cc Cu04 [LSCO] 
shows a strong growth of a contribution linear in temperature as the doping is reduced in the over- 
doped region toward optimal. This linear term is a signature of non-Fermi liquid behavior. We 
find that the appearance of a linear term in the resistivity can arise in a semiclassical Boltzmann 
transport theory which uses renormalized quasiparticle scattering rates and an empirical band struc- 
ture fitted to ARPES data on LSCO. The linearized Boltzmann equation is solved numerically by 
discretizing the Brillouin zone in a way that fits best to the Fermi surface geometry. The main 
trends in the development of the anomalous temperature dependence are well reproduced. There 
is a substantial underestimation of the magnitude of the resistivity which is expected in view of 
the moderate to weak values we chose for the onsite repulsion to stay within the one-loop renor- 
malization group approximation. The analysis was extended to the Seebeck coefficient with similar 
agreement with the main trends in the data. 



I. INTRODUCTION 



The unconventional temperature scaling of the in plane 
resistivity in high-T c cuprates has been of considerable 
interest for many years. It is now widely accepted that 
overdoped cuprates exhibit a T-linear resistivity at high 
temperature which crosses over to a T 2 dependence at 
lower temperature [IHl]. The onset of the linear tem- 
perature dependence appears to coincide with the onset 
of superconductivity and the crossover temperature to 
a conventional T 2 dependence drops monotonically with 
the doping level to zero, roughly around optimal doping. 
Hussey and collaborators measured the resistivity on a 
number of cuprates particularly the single layer T12201 
in the overdoped region and LSCO into the underdoped 
region. Very recently the Barisic et al [4] reported results 
for the resistivity for more cuprates, e.g. Hgl201 and 
YBCO over a wider density range. They found that when 
the resistivity was normalized to a give a resistance per 
CU4O4 plaquette, pn, universal results for all cuprates 
followed with identical coefficients. Further they argued 
for a universal trend in both the linear (p\j « A^T) and 
quadratic (/?□ ~ A 2 uT 2 ) regimes, with the coefficients 
A\ and A2 proportional to the inverse hole density. 
As Barisic et al [4] remarked, if one uses a Drude formula 
for resistance, this behavior is indicative of a carrier den- 
sity that is proportional to the density of holes rather 
than electrons. Such behavior is consistent with a doped 
Mott insulator scenario. Anderson [5HZ] has argued that 
the cuprates as doped Mott insulators should be in the 
strong coupling regime of the Hubbard model at all hole 
densities. He has put forward the Hidden Fermi Liq- 
uid [HFL] ansatz to describe the connection between the 
Fermi surface before Gutzwiller projection and the elec- 
tronic state that results from the removal of doubly oc- 
cupied states by Gutzwiller projection. This HFL ansatz 
gives a resistivity with the observed temperature depen- 
dence, crossing over from a linear to a quadratic form as 
the temperature is lowered. 



An earlier viewpoint interprets the cuprate phase di- 
agram in terms of an underlying quantum critical point 
separating the overdoped and underdoped regions. In 
these theories the effective interactions between electrons 
are mediated and controlled by the order parameter fluc- 
tuations associated with the quantum critical point [8 
and [9]. Another approach recently put forward by Kokalj 
and coworkers [10] is based on an anisotropic marginal 
Fermi liquid self energy. Their phenomenological theory 
gives good fits to the resistivity and also to the angular 
dependent magnetoresistance [ADMR] results [11] and a 
number of other properties. 

In an earlier study, which employed a functional 
renormalization group [FRG] treatment of the two- 
dimensional Hubbard model to describe the overdoped 
region, Ossadnik et al [12 found good qualitative agree- 
ment with the experimental results on overdoped T12201 
for the temperature and angular dependence of the relax- 
ation rate determined from the magnetoresistance. Since 
overdoped T12201 shows a conventional full Fermi sur- 
face in ARPES at zero magnetic field [13 j and in quan- 
tum oscillation experiments at finite magnetic fields [14] 
a perturbative treatment is of interest. Note that the 
FRG calculations are evaluated only to one-loop order 
and therefore are quantitatively reliable only for weak to 
moderate values of the Hubbard onsite interaction, U. As 
a result, a quantitative comparison to the experimental 
data shows a substantial underestimate of the magnitude 
of the calculated temperature dependent resistance, e.g. 
see the comparison in Kokalj and McKenzie [T5] . 

In their calculation of unusual transport behavior, Os- 
sadnik and collaborators analyzed the quasiparticle scat- 
tering vertex for overdoped cuprates within a functional 
renormalization group (FRG) approach [12]. They found 
an anisotropic term in the quasiparticle scattering rate 
with an angular form similar to the experiment and with 
an unconventional T-linear dependence which increased 
when approaching optimal doping starting from the over- 
doped side. They argued that this unusual term appear- 
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FIG. 1. Experimental results of Cooper et al. [2]. The in-plane resistivity of LSCO is plotted as a function of the temperature 
for samples with doping levels from p = 0.17top = 0.33. The red line represents data without external magnetic field, cf. the 
superconducting transition, while the red diamonds represent data in high magnetic field which suppresses the superconductivity. 
The blue dotted line is a fit with to the normal state resistivity extrapolated down to zero temperature. 



ing in the imaginary part in the quasiparticle self-energy 
could lead also to an equivalent temperature dependence 
for the in-plane resistivity. 

The goal of this study is to expand this idea and move 
beyond a self energy analysis of the transport proper- 
ties to a full calculation of transport properties using 
solutions of the Boltzmann equation. This is possible 
within a FRG approach since the method computes the 
strongly renormalized quasiparticle interactions. We de- 
termine the matrix elements for the scattering events 
within the FRG framework presented in [12 and use 
numerical techniques to solve the linearized Boltzmann 
equation. A detailed numerical treatment is required in 
view of the strong anisotropy in the transport life times 
due to both the anisotropic quasiparticle interactions and 
the strongly angle dependent Fermi velocities introduced 
by the anisotropy of the band structure. The importance 
of full transport calculations was stressed by Hlubina and 
Rice in an early paper [T6] . 

As will be discussed later, our calculations agree qual- 
itatively with the unconventional temperature depen- 
dence but the magnitude of the resistance is substantially 
lower. This disagreement in the absolute values is in line 
with the restrictions on U required to justify the one-loop 
approximation in the FRG. Another discrepancy between 
our numerical results and the experiment is the lack of 
scaling of the resistance coefficients, Aiu and A 2 n, with 
1/p the inverse hole density as reported by Barisic et al 
[4]. Note however, our results are for the overdoped sam- 



ples. The data in Fig. 5 of Ref. [4 deviate from the 
inverse doping scaling in the overdoping density range. 
In the case of the linear coefficient, A^u the data show 
a linear decrease and vanish at the onset of supercon- 
ductivity in line with the earlier results of Hussey and 
coworkers [2]. 

A good example of such an onset of a non- Fermi- liquid- 
like T-linear contribution to the resistivity can be found 
in the resistivity measurements by Hussey and cowork- 
ers [2] upon reducing the hole concentration from the 
overdoped regime toward optimal doping in LSCO (cf. 
Fig. [I]). They observed a strong growth of the tem- 
perature dependent contribution to the resistivity with a 
decrease in the hole doping together with a remarkable 
trend toward linear resistivity that goes down to lowest 
temperatures. They used strong magnetic fields to ex- 
plore the normal state transport properties within the 
superconducting dome. The analysis of the temperature 
dependence of the resistivity was based on a second-order 
polynomial fit to the experimental data over a wide tem- 
perature range: p(T) = -f aiT + o^T 2 . This allows to 
distinguish doping regimes with dominant linear versus 
dominant quadratic temperature dependence. 

We also solve the Boltzmann equation to obtain the 
Seebeck coefficient which describes the charge transport 
in response to a thermal gradient. Standard Fermi liq- 
uid theory predicts a metallic Seebeck coefficient with 
a linear T-dependence. The experiments by Laliberte 
and coworkers [17] showed substantial deviations from 
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the standard form with a complex density dependence. 
Again our calculations reproduce the qualitative trends 
in the temperature and density dependence very well. 

In Sec. [IT] we show our numerical results and compare 
with the experimental data of Refs. [2H4]. We also ex- 
amine the validity of Matthiessen's rule which relies on 
a decoupling of impurity and two-particle scattering in 



Sec. |III| Finally in Sec. [TV} we study the thermoelectric 
effect of LSCO, use the Seebeck coefficient as a sign of 
non Fermi liquid transport and examine the evidence for 
quantum critical behavior. 

The numerical method that we have developed to 
study the transport properties of correlated materials 
such as LSCO is explained in detail in the appendix, 



App. |A]|Cl We solve the semiclassical Boltzmann equa- 
tion taking the full angular and energy dependence of 
the distribution function into account. The two impor- 
tant inputs of the model are the explicit form of the band 
structure and the temperature dependent calculation of 
the quasiparticle scattering rates. 



II. NUMERICAL RESULTS 

We begin by reviewing our numerical results. We use 
a semiclassical model to describe the normal state charge 
transport of overdoped LSCO. For the quasiparticle dis- 
persion we used a two-dimensional tight binding model 
and a doping dependent hopping parameter renormaliza- 
tion scheme based on the values determined in Ref. [18] 
from the Fermi surface in APRES data on LSCO. The 
use of phenomenological input for the quasiparticle dis- 
persion is important because the renormalization of the 
hopping integrals is of two-loop order in the FRG anal- 
ysis that we use to determine the scattering vertex. The 
calculation of two-loop diagrams for the RG flow is, how- 
ever, beyond the scope of this study. For details about 
the dispersion, cf. App. [Aj The FRG technique that we 
use to compute the renormalized quasiparticle scattering 
rates is outlined in App. [B] together with brief overview 
of the doping and temperature dependence of the scatter- 
ing vertex and a discussion of the parameters that enter 
the RG calculation. 

Taking the quasiparticle dispersion and the scattering 
vertex as input, the conductivity is computed by solving 
the linearized Boltzmann equation with its full angular 
and energy dependence. Due to the strong anisotropy 
in the scattering rates and the quasiparticle velocities, 
a full solution of the Boltzmann equation rather than a 
single relaxation time approximation is required. The 
collision integral is computed by introducing an efficient 
discretization of the Brillouin zone and the Boltzmann 
equation is then solved numerically. This method of solv- 
ing the Boltzmann equation and discretizing the Brillouin 
zone is explained in App. [C] 



A. Simulation estimates of p(T) for LSCO in a 
range of doping 

The in-plane resistivity in our model of the normal 
state of overdoped LSCO is shown in Fig. [2] (data for 
U\ = 2.0 1). We also plot the resistivity for different 
initial values for the unrenormalized onsite interaction 
parameter U. The data corresponding to the different 
values of U\ are derived from renormalized scattering 
vertices while the data for [Too are computed for particles 
interacting via the bare onsite repulsion. For a compar- 
ison of these calculations, cf. Sec. II C The values for 
the hole-doping p are chosen to match the experimental 
measurements of Fig. [I] 

Our results compare well to many qualitative aspects 
of the experimental data of Ref. [2 . First, there is the 
strong growth of the resistivity when reducing the hole 
doping from strongly overdoped values of p ~ 0.35 to 
optimal doping around p ~ 0.17. Second, there is clear 
evidence for a growing linear term in the temperature 
dependence of p around optimal doping. For large val- 
ues of the hole doping, p(T) displays a rather quadratic 
temperature dependence. On the other hand, for dop- 
ing concentrations below p « 0.25 a linear T-dependence 
becomes obvious. We analyze our data in the same 
scheme as in the experimental paper [2 using f its with 
the form p(T) = a + a x T + a 2 T 2 (Sec. |IIB1[ ) or also 
by introducing a "local scaling exponent x(T,p), 
p(T) = a + aoo^ (Sec. |HB~3). 



e.g. 



For a quantitative comparison to the experiment we 
convert our two-dimensional data to the experimental 
units of a three-dimensional sample (cf. App. |P| ). From 
Fig. [2] we find that our calculated values for p(T) at 
T « 200 K are roughly by a factor four or five smaller 
than in experiment (cf. Fig. [TJ. An underestimation 
of the resistivity is consistent with the fact that our un- 
renormalized interaction strength U = 2.0 1 is consider- 
ably lower than expected for LSCO. 

In spite of the weak coupling approximation in our cal- 
culation of the scattering rates, the resistivity is in good 
qualitative agreement with the observed temperature de- 
pendence. 



B. Scaling behavior 

1. T + T 2 fit to the resistivity curves 

Next we analyze the numerical resistivity data using 
the polynomial fit up to second order p(T) = -\-a\T-\- 
a 2 T 2 1 in analogy to Cooper et al, cf. Ref. [2 . They fit 
their resistivity data over the entire temperature range 
and obtain the fit coefficients ql\ and a 2 as a function 
of the hole doping, cf. Fig. [3] (left panel). They found 
the linear coefficient ot\ grows strongly with a decrease 
in doping while there is a constant or decreasing value of 

OL 2 . 
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FIG. 2. Numerical calculation of the normal state resistivity of LSCO (black line, Ua = 2.0 1) as a function of the temperature 
for several doping levels: The resistivity is plotted for several values for the initial onsite repulsion U . The RG cutoff is given 
by A = 0.01 1, except for the dashed line (Uoo) computed from unrenormalized scattering rates. The T — » extrapolation of p 
corresponds to the residual resistivity due to impurity scattering (ni m pW|_ im p = 0.25 1 2 ). We have converted the temperature 
and resistivity to the units used in the experiment. For details cf. App. [D] 



The linear ol\ and quadratic ol^ coefficients from fit- 
ting our numerical results are shown in Fig. [3] (right 
panel) for the entire range of hole doping from p « 0.15 
to p « 0.35. The fit was performed for temperatures 
ranging from 30 K up to 300 K. For the comparison of 
the two coefficients we introduce an average temperature 
scale T m = 150 K in the conversion to resistivity units. 
The coefficient representing the residual resistivity, ao, is 
discussed later. 

The direct comparison of our fit coefficients to the ex- 
periment yields very good qualitative agreement. The 
linear term shows a strong increase with reduced hole 
concentration toward optimal doping while above p « 
0.30 it becomes indistinguishable from zero. On the other 
hand, the quadratic contribution does not grow and even 
decreases almost down to zero toward optimal doping for 
larger U. A vanishing ot2 indicates a within our model a 
purely linear temperature dependence. Whether or not 
the quadratic coefficient of the fit to the experimental 
data vanishes is not really clear from Ref. [2 j, although 
it appears that the ot2 determined with the conventional 
T + T 2 fit indeed decreases and finally vanishes around 
p & 0.19, while the quadratic coefficient determined from 
a parallel resistor fit remains constant. 

The coefficients in the overdoped regime disagree with 
the reported 1/p dependence by Barisic et al. [4 . As 
already mentioned in the introduction, in the overdoped 
regime, their data suggest a deviation from the remark- 
able 1/p dependence of the linear term in the underdoped 
regime. The experimental data of both Ref. [2] and [4] 
indicate that the linear term eventually goes to zero for 
sufficiently large hole density, p ~ 0.30. 



2. Scaling regimes 

An alternative analysis of the experiments was re- 
cently put forward by Hussey et al. [3 j. Taking 
the derivative p' = dp(T)/dT they distinguished two 
regimes which are not immediately obvious in the re- 
sistivity data, but in which p(T) scales differently with 
temperature. Hussey and coworkers identified a high- 
temperature regime where p' is constant, i.e. the resistiv- 
ity depends linearly on T. Toward lower temperatures p' 
deviates downward from this constant value. With this 
rather abrupt change one can associate a temperature 
T co h, called "coherence temperature" in Ref. [3 which 
decreases with lowering the hole concentration and dis- 
appears probably around the same value of p where the 
characteristic temperature for the pseudogap vanishes ex- 
trapolated from the underdoped side. In Fig. [4] the dop- 
ing dependence of T co h measured in Ref. [3] is shown. 

Plotting p'(T) for different doping levels in our cal- 
culation, we can analogously identify distinct regimes in 
the temperature scaling of the resistivity, cf. Fig. [5] (left 
panel). There is a low-temperature regime with finite 
slope of p'(T), an intermediate temperature regime with 
finite, but smaller slope, and a high-temperature regime 
where p'(T) is constant. The temperatures that separate 
these three regimes will be referred to as T\ and T2. 

The temperature scales T\ and T2 are displayed in Fig. 
[6] The doping dependence of T2 shows a very similar 
trend as Tqq^. This is interesting because there is no 
real loss of quasiparticle coherence in our semiclassical 
model based on the Boltzmann transport equation. Note 
that we have not considered the values for T\ and T2 for 
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FIG. 3. left: Experimental results obtained by Cooper et al. 2 . The quantities a\ and «2 are coefficients of fit to the 
experimentally measured resistivity (Fig. [TJ: black filled squares: fit as in our numerical analysis, p — ao + ol\T + a2T 2 -fit; 
green filled circles: parallel resistor fit, 1/p = l/(ao + a\T + CX2T 2 ) + l/p m ax; open symbols correspond to the same analysis 
for experimental data of Ando et al. ^9] right: Linear (a±) and quadratic (0.2) term of the p — ao + ol\T + (X2T 2 fit to our 
numerical resistivity data of Fig. [2] as a function of the hole-concentration p: Shown are the coefficients for different initial 
values of the onsite repulsion U and for the calculation with unrenormalized interactions (Uoo)- We convert the coefficients to 
resistivity units by introducing the average temperature scale T m = 150 K of the fit regime. 



p < 0.2 as it becomes difficult to distinguish the regimes 
for p^ 0.19 

It is interesting to relate the onset of the unconven- 
tional scaling, T\ , to the proximity of the Fermi energy to 
the van Hove level, £? v h- We introduce the temperature 
scale Ta vH = 4|/i — -EvhI/^b? which is shown as the red 
squares in Fig. [6] The factor of 4 is introduced for con- 
venience, as we consider only those low-energy states rel- 
evant for the transport, that live within an energy range 
of 4:ksT from the Fermi level. The doping dependence of 
Ta vH and 7\ is very similar and we conclude that prox- 
imity to a van Hove singularity plays an important role 
in the development of non- Fermi liquid behavior. 

To compare to the analysis in Ref. we fit the high- 
temperature linear part of p' with a constant ai(oo) 
and the low-temperature regime by a linear function 
ai(0) + 0L2F . In the high-temperature regime, ol\[oo) 
represents the coefficient of the dominating linear term 
in the resistivity, while ai(0) corresponds to the coeffi- 
cient of the dominating linear term in the limit T — >• 0. 
The results for these coefficients are shown in Fig. [5] 
(right panel). 



3. Derivation of the local scaling exponent 

A different way to illustrate the regime change on the 
resistivity data is to define an exponent of p(T) for each 
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FIG. 4. Experimentally determined coherence temperature 
Tcoh (p > 0.19) as a function of the hole doping of LSCO. The 
coherence temperature is where the temperature derivative 
of the conductivity deviates from its constant value at high 
temperature (data adopted from Ref. [3])- 



temperature formulated as 

p(T) = a + aT x , 



x = x(T). 



(1) 



The exponent x(T) is obtained by taking the logarithmic 
derivative of the p(T), 



x(T) = 



d 



d\uT 



ln(p(T) - a ), 



(2) 
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FIG. 5. left: Plot of the temperature derivative of the resistivity p'(T): Three intervals are indentified in which p(T) shows 
different scaling behavior. Shown is p' (T) for four different values of the hole doping illustrating the doping dependence of T\ 
and T 2 . For the bottom right inset at doping level p « 0.20, the temperature scales cannot be clearly identified any more. 
In Fig. [6] Ti and T 2 are plotted as a function of the doping, right: Limiting behavior of the linear term that dominates 
the temperature dependence of the resistivity: The coefficient ai(0) represents the T — » extrapolation of the slope in the 
resistivity curve while ai(oo) corresponds to the coefficient of the dominating linear term in the high-temperature regime. The 
strong growth of ai(0) for doping levels p < 0.25 reflects the fact that the unconventional scaling behavior extends down to low 
temperature for these dopings. The definition of oa(0, oo) only allows these quantities to be obtained in a finite doping range. 



where ao denotes the residual resistivity. The results for 
x(T) over a wide range in doping are shown in Fig. [6] 

Even though such a local scaling law has large uncer- 
tainties, it is intriguing to find a very good qualitative 
match between the function of x(T,p) and the generic 
scaling behavior of the cuprates, cf. and Ref. pQ. In the 
overdoped low-temperature region of the phase diagram, 
the local exponent rises to the value of x « 2 correspond- 
ing to a standard Fermi liquid. Around optimal doping, 
x(T,p) is always well below 2 down to the lowest temper- 
ature of our study. In the high-temperature regime, the 
onset of unconventional scaling (x(T,p) < 2) is shifted to 
higher temperature with increasing doping. 



ar resolved scaling analysis 



life times. The scattering rates that cause these short 
life times are strongly enhanced as the doping level de- 
creases toward optimal doping. To examine the role of 
strong scattering of the anti-nodal quasiparticles for the 
unconventional temperature dependence of the overall re- 
sistivity, we interpret the angular patches of our Brillouin 
zone discretization (cf. Fig. 13 in the appendix) as in- 



dividual resistors and parametrize the low-energy quasi- 
particle states as a parallel resistor network. Each resis- 
tor is characterized by the conductivity a* and resistivity 
pi = l/<Ji where the index i labels the the angular patch. 
The total conductivity and resistivity are determined by 



a 



= Eva 



(3) 



Using the non-equilibrium distribution function, we 
can analyze the angular dependence of this unconven- 
tional scaling and compare to the difference in the tem- 
perature dependence between isotropic and anisotropic 
scattering rates in the FRG study [T2] . The quasiparti- 
cle scattering rates are strongly anisotropic and especially 
the quasiparticles close to the saddle points (anti-nodal- 
direction) in the Brillouin zone have very short transport 



In Fig. [7|we depict the the angle-resolved scaling by plot- 
ting the linear and quadratic contributions to the patch 
resistivity as functions of the angular patch index i and 
the hole concentration. The temperature dependence of 
each resistor was analyzed by a linear plus quadratic fit, 
pi(T) = Oq + a\T + a l 2 T 2 . In order to compare linear 
and quadratic term, we give them a weight factor that 
represents their contribution to the total resistivity and 
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FIG. 6. Temperature and doping dependence of the scaling 
exponent x(T,p) denned in Eqs. |l]) and |2| for our model 
of LSCO in color-code. The temperature scales T\ and T2, 
cf. Sec. |IIB 2| are shown as black squares. The temperature 
scale Ta vH corresponds to the energy difference between the 
chemical potential and the energy level of the van Hove sin- 
gularity, /cbTa vH = 4\/n — E v ii\ (we consider low-energy states 
up to an energy of 4k bT around the Fermi energy). Note 
that the unconventional scaling x < 2 goes down to the low- 
est temperatures at optimal doping. The Lifshitz transition 
is indicated as the horizontal dashed line. 
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FIG. 7. Comparison of linear vs. quadratic temperature scal- 
ing of the resistivity for the angular patches (cf. Fig. 13) 



On top, we have plotted the weighted linear term A\, and the 
lower plot shows the weighted quadratic term A|, defined in 
Eq. (pi. These A^s are weighted according to the contribu- 
tion of the angular patch i to the total charge transport. The 
patches range in angular direction from = (patch 1) to 
= 7r/2 (patch 10). The A l a 's are evaluated at the tempera- 
ture Tm — 150 K, the average temperature of the fit region. 



normalize their sum, the total resistivity p, to unity, 
p al + a\T + a\T 2 



with 



Ai 



Pt 



al T a 



STAi + A\ + Ai (4) 



0,1,2. 



(5) 



Thus A l a represents the normalized contribution of the 
angular patch i to the term in the total resistivity that 
scales with temperature exponent a. 

From Fig. [7] one can clearly see that in the strongly 
overdoped regime only the quadratic term is present and 
is essentially independent of the angle. On the other 
hand, below p w 0.30, the linear term grows showing 
strong anisotropy. The linear term appears first close 
to the nodal direction (patches 4 to 7), and constantly 
spreads out. Below p ~ 0.27 the linear term starts to 
dominate the transport. The quadratic term, however, 
remains finite in the nodal direction down to the lowest 
doping levels. 

It is interesting to see that around optimal doping, 
where the Fermi surface is very close to the van Hove sin- 
gularity, the charge transport from the angular patches 
1 and 10, that contain the saddle points, essentially van- 
ishes. This is consistent with the strong growth of the 



scattering rates from the FRG analysis. These angu- 
lar patches are short-circuited as a consequence of strong 
scattering, cf. Ref. [16] . Our discussion here is fully com- 
patible with the behavior of the angle-dependent scat- 
tering rates found by the FRG study by Ossadnik and 
collaborators [12]. 

We can also analyze the current distribution as a func- 
tion of the temperature. In Fig. [8] we display the angular 
distribution of the charge current for two different dop- 
ing levels. At the lowest temperatures the distribution 
function changes its angular dependence strongly and the 
weight of the current is shifted from small angles to larger 
angles, while above a certain temperature the angular 
distribution remains invariant. 

These regimes can be associated with two different 
kinds of scattering mechanisms. While at low tem- 
perature impurity scattering dominates, for sufficiently 
high temperature the transport is mainly determined by 
particle-particle scattering. Since only umklapp scat- 
tering contributes noticeably to the resistivity, the two- 
particle scattering contribution is anisotropic, while the 
^-potential impurities give isotropic scattering. This dif- 
ference accounts for a qualitative change in the angular 
dependence of the distribution function. We will return 
to the interplay of the two scattering mechanisms in the 



context of a discussion of Matthiessen's rule, Sec. Ill 



8 




FIG. 8. Plot of the angular distribution of the current for doping levels p = 0.23 and p — 0.29. The curves correspond to 
temperatures between 30 and 600 K indicated by their color. Shown is merely the angular range between and tt/2 due to the 
four-fold lattice symmetry. 



C. Effect of the vertex renormalization 

We turn now to the role of the renormalization of 
the scattering vertex in the unconventional nature of the 
transport properties. 

As explained in App. |B 2| we used a value of U\ = 2t 
for LSCO. When we extend our calculation to a series of 
lower initial values of U\ the temperature dependence of 
the resistivity and the coefficients of the corresponding 
linear plus quadratic fit are shown Fig. [2] and in Fig. [3] 
In addition we show the resistivity for the unrenormalized 
interaction U^ = 2.0 1 which corresponds to an FRG 
treatment with infinite cutoff. In this case the scattering 
rates are determined only by the scattering phase space 
which is specified by the Fermi function and the geometry 
of the Fermi surface. 

Comparing the resistivity curves for U\ = 2.0 1 and 
Uoo = 2.0 1, we observe very similar overall trends as 
function of doping. The high-temperature resistivity 
grows upon reducing the hole-concentration and the T- 
linear part increases at the same time. It is obvious 
that the scattering phase space plays an important role 
in the change of the behavior of resistivity for different 
hole-concentrations. Overall the behavior is, however, 
less pronounced for the unrenormalized than renormal- 
ized vertex as is clear in Fig. [3] The renormalization of 
the vertex causes a sizable increase of the resistivity as 
one can see in the comparison of the curves for U\ > 2.0 1 
and Uoo — 2.0 1, and it also strongly enhances the trend 
toward linear resistivity. This feature is also easily iden- 
tified from the [/-dependence of the quadratic term in 
Fig. [3j For stronger initial interaction, and therefore 
stronger renormalization effects, the quadratic term is 
substantially reduced, even as the total resistivity, the 
sum of linear and quadratic term is increased. Thus, it 
is evident that correlation effects included in the vertex 
renormalization by the FRG scheme, make a substantial 
contribution of the unconventional transport properties, 



in particular, in the regime close to optimal doping. 

The vertex renormalization leads to a non-trivial mo- 
mentum space structure of the effective scattering vertex 
which can be translated to a longer ranged real space in- 
teraction compared to the local Hubbard-/7. As a second 
effect of the FRG analysis, the momentum space struc- 
ture but also the absolute values of the scattering vertex 
become temperature dependent. In order to discriminate 
between these two effects we modified our calculations in 
a way as to include only the momentum space renormal- 
ization of the vertex, but neglect the temperature depen- 
dence. This was implemented by stopping the RG flow 
with a large cutoff A. In this case large refers to the cut- 
off exceeding the largest temperature of our study. While 
this completely suppresses the temperature dependence 
of the scattering vertex, it still introduces momentum 
space anisotropy. 

The result resembles the behavior of the simple un- 
renormalized onsite repulsion U^ = 2.0dt of Fig. [2] The 
introduction of momentum anisotropy through the RG 
flow has almost no qualitative effect on the temperature 
dependence of the resistivity. From this comparison we 
can conclude that the pronounced linear resistivity in 
the limit of T — >• is a consequence of the strongly tem- 
perature dependent quasiparticle scattering rates. This 
supports the idea of Ref. [12] where the unconventional 
temperature dependence of the scattering rates through 
the RG flow were described. 

At high temperature, however, the growth of the 
linear-temperature term when the doping is decreased 
toward optimal doping, is mainly due to the Fermi sur- 
face geometry and the temperature dependence of the 
scattering volume, while the momentum anisotropy of 
the scattering interaction is responsible for the enhanced 
value of the resistivity. 
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III. MOMENTUM RELAXATION MECHANISM 
AND MATTHIESSEN'S RULE 

In the cuprates the resistivity is caused by electron- 
electron interactions which raises the question, how the 
total momentum of the entire distribution function is re- 
laxed. In the following we will discuss which contribu- 
tions are combined in the resistivity calculations of the 
previous section. 



A. Impurity and umklapp scattering 

In impurity scattering processes, the momentum dif- 
ference between initial and final state is transferred to 
the lattice and determines the residual resistivity as the 
zero-temperature limit. A deeper analysis shows that the 
resistivity from impurity scattering can have a weak tem- 
perature dependence due to thermal broadening of the 
distribution function and a possible bias toward higher 
or lower quasiparticle velocity. In our calculation the im- 
purities are modeled as ^-potentials. 

A second and more complex contribution to the electri- 
cal resistivity originates from electron-electron scattering 
with finite momentum transfer to the lattice, so-called 
umklapp scattering, where a momentum corresponding 
to a reciprocal lattice vector is transferred to the lattice. 
Note that the Fermi surface has to be sufficiently large 
in order to allow for umklapp assisted scattering between 
different momenta on the Fermi surface. 

For a square lattice we can distinguish between two 
kinds of umklapp scatterings: with momentum transfers 
of the type (27r/a, 0) and (27r/a, 2ir/a) to the lattice. The 
first kind of umklapp scattering has a large phase space 
for a sufficiently large Fermi surface with regions where 
the x-component of the Fermi vector satisfies the condi- 
tion \kp x \ < K / a and other regions with \kp x \ > tt/cl. 
The second type of umklapp process is only possible, if 
the Fermi surface intersects with the so-called umklapp 
surface, which is the rotated square (diamond) that con- 
nects the four saddle points at (0, d=7r/a) and (±7r/a,0). 
This umklapp scattering has relatively small phase space 
due to the stringent momentum constraints. Within our 
model for LSCO the (27r/a,0) umklapp process largely 
dominates over the (27r/a, 2it/a) umklapp process for all 
doping levels as far as the charge transport resistance is 
concerned. On the strong overdoped side, (27r/a, 2ir/a) 
umklapp scattering is even completely absent in the low- 
temperature limit. 

Due to the geometrical constraints for umklapp scat- 
tering the momentum space structure of the scattering 
vertex can play an important role. The temperature de- 
pendence of the resistivity is controlled by umklapp scat- 
tering. If the scattering rates become especially high on 
certain parts of the Fermi surface, their contributions to 
transport will be suppressed yielding the possibility for 
unconventional temperature scaling. 



B. Deviations from Matthiessen's rule 

Often different contributions to the resistivity are as- 
sumed to be decoupled and one adds them as serial resis- 
tors. This decoupling leads to Matthiessen's rule, which 
is only justified, if each individual contribution can be 
well described by a single relaxation time. In this case 
the collision integral in the Boltzmann equation directly 
splits into the scattering channels. In the general case, 
when such a single relaxation time approximation is in- 
sufficient, the collision integral does not decouple. 

In Fig. [8] we have seen that in the low-temperature 
regime the angular distribution of the current changes 
its form. This was attributed to a crossover from an 
impurity scattering dominated regime to a two-particle 
scattering regime. These different scattering mechanisms 
lead to different current distributions due to their differ- 
ence in the scattering geometry. While impurity scat- 
tering is isotropic, the umklapp processes responsible for 
the momentum relaxation in the two-particle scattering 
channel are highly anisotropic leading to a non-trivial 
interplay of impurity scattering with two-particle inter- 
actions and a potential breakdown of Matthiessen's rule. 

To examine the deviations from Matthiessen rule we 
compare our calculation with the sum of the "individ- 
ual" contributions to the resistivity. First we remove 
the impurities to obtain a "bare" electron-electron resis- 
tivity pee, second we estimate pi mp by considering free 
(non-interacting) electrons scattered by impurities. 

In a first step we compare the linear and the quadratic 
coefficients of a T + T 2 -fit for the clean limit, p ee , to 
the coefficients in the original calculation including im- 
purities of Fig. [3j The quadratic coefficient grows weakly 
down to p « 0.25 and remains constant at p < 0.25, while 
for the full calculation the quadratic coefficient decreased 
below p ~ 0.30 and even went to zero for sufficiently 
strong coupling. Furthermore, the linear term is some- 
what smaller without impurities. This difference in the 
behavior of the coefficients is a violation of Matthiessen's 
rule. 

In order to make this more evident, we plot the devi- 
ation from Matthiessen rule, 5m, quantitatively defined 
as 

$M = [Pfull - (Pee + Pimp)] Mull j (6) 

as a function of the temperature for different doping lev- 
els. The result is shown in a color plot in Fig. [9] We im- 
mediately observe that the deviation from Matthiessen 
rule is strongest at the optimally doped regime at low 
temperature. The difference in the scattering geome- 
try for impurity scattering and two-particle scattering 
is most pronounced in this part of the diagram. Here the 
deviation 5m reaches values as large as 16%. 

It should be noted that the one-loop FRG approxi- 
mation is less reliable at low temperature and low dop- 
ing, as pointed out earlier. Hence the magnitude of 5m 
in our calculation carries considerable uncertainties in 
this part of the diagram. Nevertheless, the violation of 
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FIG. 9. Deviation from Matthiessen rule: We subtract the 
sum of the resistivity for pure impurity scattering and pure 
electron-electron scattering from the resistivity of the full cal- 
culation pf u n and normalize it to pfuii- It is clearly visible 
that the deviation 5m (cf. Eq. |6])) is most pronounced in the 
critical region of the phase diagram. 



Matthiessen's rule is an generic property for systems fea- 
turing a crossover from isotropic impurity scattering to 
strongly anisotropic inter-particle scattering dominated 
transport. 



IV. THE SEEBECK EFFECT IN OVERDOPED 
LSCO 

The Seebeck effect is a further example which allows 
us to test Fermi liquid properties in transport. This ef- 
fect as the electric response to a temperature gradient is 
based on thermo-diffusion of charge carriers and can be 
analyzed with Boltzmann transport theory too. 



A. Experimental observation of the Seebeck effect 
in Eu-LSCO 



Laliberte et al. measured the Seebeck coefficient Q of 
La1.8_sEuo.2Sra.CuO4 [Eu-LSCO], cf. Ref. [17]. While 
their study aims at probing the Fermi surface reconstruc- 
tion scenarios in the underdoped regime of the phase di- 
agram, they also present data up to a hole-doping level 
of p = 0.24. Data of the Seebeck coefficient are displayed 
in Fig. [To] (left panel) for the doping regime p > 0.11. 
Europium doping of LSCO is believed to have a stronger 
influence on the properties of underdoped samples, but 
little effect is expected in the overdoped regime. 

For underdoped samples Q/T shows a sign change, 
starting negative at low temperature it turns positive at 



higher temperature. On the other hand, for overdoped 
samples Q/T remains positive for all temperatures mea- 
sured and shows a rather strong increase toward low tem- 
perature near optimal doping. The inset of Fig 
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(left 

panel) demonstrates the log(l/T) dependence of Q/T at 
low temperature, which has been interpreted as a signa- 
ture of quantum criticality. 



B. Numerical results 



We have numerically calculated the Seebeck coefficient 
for LSCO within our semiclassical approach, as explained 
in App. 



C 2 The results are shown in Figs. 10 (right 



panel) and 11 A good overview on the behavior of the 
Seebeck coefficient as a function of doping and tempera- 
ture is obtained from Fig. 
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(left panel). At the doping 
p c « 0.19 we observe a sign change as we cross the Lif- 
shitz transition at low temperature. Interestingly, for 
p < p c increasing the temperature leads to a sign change 
from negative at low to positive at higher temperature. 
Turning to the other limit of strong overdoping the See- 
beck coefficient turns negative again as expected for elec- 
trons is a band with small filling. We also observe a 
doping dependent temperature scale where Q displays a 
rapid increase. Interestingly, this temperature compares 
well with T\. 

In order to compare our numerical with the experi- 
mental results, we have plotted Q/T in Fig. 10 for some 
selected values of hole doping. The inset shows the low- 
temperature regime for the doping level p = 0.21 on a 
logarithmic temperature scale. 

The general structure of the calculated Seebeck coef- 
ficient is in good qualitative agreement with the exper- 
imental data measured for Eu-LSCO Fig. 10 In both 



cases, there is a critical doping level below which the 
Seebeck coefficient changes sign with increasing temper- 
ature, starting from negative in the T — >• limit. Above 
this critical doping level, Q/T becomes very large for 
low temperature and strongly decreases with tempera- 
ture. The low-temperature rise of Q/T becomes less 
pronounced with increasing hole doping. In the low- 
temperature regime, Q/T acquires a log(l/T) depen- 
dence which is clearly visible in both the experimental 
data and the numerical data below 60 K. 

A small discrepancy between our model and the exper- 
iment is the value of the critical doping level. We have de- 
signed the band structure that we use to describe LSCO 
such that the Lifschitz transition is located at p c « 0.19, 
cf. App. [A] The van Hove singularity at the Lifschitz 
transition is the origin for the critical behavior in our 
model and thus constrains the critical doping level to 
p c . In the experimental study, however, a critical doping 
level of 0.24 is reported. 
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FIG. 10. left: Experimental data of the Seebeck coefficient Q normalized to the temperature for Eu-LSCO 
(Lai.s-^Euo^Sr^CuCU), adapted from Ref. [T7]. right: Numerical calculation of the Seebeck coefficient Q normalized to 
the temperature for our model of LSCO. Note that the temperature range does not extend to zero. 
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FIG. 11. left: Seebeck coefficient as a function of hole doping and temperature, right: Derivative of Q/T with respect to 
the temperature: Fermi liquid theory predicts a linear temperature dependence of the Seebeck coefficient. Thus, any non-zero 
value of this derivative indicates deviations from Fermi liquid theory. 



C. Signature of critical transport properties 

In the last section, we have discussed the low- 
temperature log(l/T) dependence of Q/T and the sign 
changes with doping or temperature. These properties 
of the Seebeck coefficient are often taken as indication 
for the presence of a quantum critical point [22l and [23] . 
These features are found in both the experimental data 
and the numerical simulation. We now want to locate 
the areas of critical behavior in the phase diagram. 

For conventional metals, the Seebeck coefficient is a lin- 
ear function of the temperature manifesting Fermi liquid 



behavior. Deviations from this linear dependence may be 
used as a signature for non-Fermi liquid physics. Thus, 
we analyze the the temperature derivative of Q/T which 
should vanish for Fermi liquids. Fig. [II] (right panel) 
displays this function. This value is largest at low tem- 
perature around optimal doping and a sign change ap- 
pears at the Lifshitz transition. This is the region where 
the strongest differences with the Fermi liquid picture 
emerge. From this finding we may identify the vicinity 
of the van Hove singularity and the presence of a Lifshitz 
transition as the origin of non-Fermi liquid behavior and 
the signatures of quantum criticality. 
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Interestingly, Fig. 11 (right panel) allows us also to see 



the temperature of the sharp increase of Q/T, cf. Fig. 



11 (left panel), appearing as a maximum of d(Q/T)/dT. 
This characteristic temperature T max separates the high- 
doping low-temperature regime coming closest to Fermi 
liquid behavior (T 2 -dependence of resistivity) from the 
lower-doping high-temperature regime which is charac- 
terized by the non-Fermi liquid T-linear resistivity, as 
discussed in Sec. Ill B 21 



V. DISCUSSION AND CONCLUSION 

We investigate the temperature and doping depen- 
dence of charge transport coefficients using a combina- 
tion of an experimentally determined band structure and 
a renormalized scattering matrix calculated in a one-loop 
renormalization group approximation. The band disper- 
sion matches the quasiparticle spectrum in ARPES for 
LSCO [18], and as such includes the renormalization of 
the dispersion with hole density due to interactions. The 
electron-electron scattering vertex develops momentum- 
space and temperature dependence as the hole density is 
reduced. We have chosen a moderate to weak value for 
the bare onsite Coulomb repulsion consistent with the 
one-loop approximation. We also refrain from including 
the pseudogap region which would require a solution of 
the strong coupling regime that results from divergences 
in the FRG flow. The Boltzmann transport equation is 
solved numerically with collision integrals for electron- 
electron and impurity scattering determined by Fermi's 
Golden rule. This semiclassical approach determines the 
quasiparticle distribution functions taking into account 
the full angular and energy dependence of the collision 
integral. 

We find good qualitative agreement with the main 
trends in the experiments on the in-plane resistivity of 
LSCO [2 . This includes the increase of the temperature 
dependent term in the resistivity upon decreasing the 
hole-concentration and the trend toward a linear temper- 
ature dependence of the resistivity p(T) around optimal 
doping. A parameterization of our results using a poly- 
nomial fit in T up to second order, inspired by the exper- 
imental analyses [2] and [4], gives consistency between the 
trends in experiment and our model. While the T-linear 
behavior of p(T) is often attributed to the existence of 
a quantum critical point around optimal doping [20] and 
[21], our discussion shows that this feature may also be 
understood as a consequence of the unusual temperature 
dependence of the renormalized scattering vertices and 
the increased scattering phase space around this doping 
regime, where a Lifshitz transition of the band structure 
leads the chemical potential very close to a van Hove 
singularity. The magnitude of the calculated resistivity 
is substantially lower than in the experiment, consistent 
with the moderate to weak interaction values and one- 
loop RG approximation 

Motivated by the experimental data analysis [3] we also 



identified crossover regimes in dp(T)/dT and find rather 
sharp changes in the temperature dependence. There 
are two characteristic temperatures T\ and where the 
latter denotes the lower bound of the T-linear regime 
of p(T) and Ti(< T 2 ) the upper bound for a purely 
T 2 -dependence (Fermi liquid regime). The intermediate 
temperature range may be considered as combination of 
both. The resulting phase diagram gives a very good ac- 
count of the trends of p(T) observed in the experiment. 
In particular, we can show that there is a clear correla- 
tion of the unconventional temperature scaling with the 
van Hove singularity, quantified through the effective dis- 
tance of the chemical potential from the van Hove level, 
k B T AvH = 4|/i - E vH \, cf. Fig. [6j 

By comparing the temperature dependence of p(T) for 
renormalized and bare scattering rates, we could demon- 
strate the importance of the scattering vertex renormal- 
ization. While in the strongly overdoped regime, the dif- 
ference in the numerical data determined from renormal- 
ized and unrenormalized scattering rates is not very pro- 
nounced, around optimal doping the phase space driven 
growth of a linear term is strongly enhanced due to renor- 
malization of the quasiparticle interactions. 

It also turns out that the impurity scattering is impor- 
tant in this respect as well, contributing to the uncon- 
ventional temperature dependence. This is surprising as 
impurities contribute usually only a weakly temperature 
dependent offset to resistivity. We found deviations from 
Matthiessen's rule in large parts of the considered phase 
diagram. This deviation is especially pronounced as we 
approach optimal doping at intermediate to low tempera- 
ture due to the different momentum dependence of impu- 
rity and electron-electron scattering. Matthiessen's rule, 
however, is recovered in the more conventional, strongly 
overdoped regime at low temperature. 

We have also studied the Seebeck coefficient within our 
model. We find qualitative agreement between our nu- 
merical data and the main trends in the measurements 
performed on a slightly different compound, Eu-LSCO. 
We believe that the Eu-doping does not strongly influ- 
ence the Seebeck coefficient in the optimal to overdoped 
doping regime which we are interested in. The calculated 
Seebeck coefficient, also displays a critical region in the 
phase diagram with pronounced deviations from Fermi 
liquid behavior. We have related the critical behavior to 
proximity of the van Hove singularity to the Fermi energy 
which is linked to the Lifschitz transition. 
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Appendix A: The quasiparticle dispersion of LSCO 

We use a band structure parametrized within a tight- 
binding model on a square lattice, including both nearest 
and next-nearest neighbor hopping 

Sk = — 2t(cos k x + cos k y ) + 4t' cos k x cos k y 

- 2t"(cos 2k x + cos 2k y ). (Al) 

The hopping parameters t, t" ', and the chemical po- 
tential /i define the geometry of the Fermi surface and 
the filling of the band. 

Based on their ARPES study, Yoshida et al. [18] have 
fitted the low-energy quasiparticle spectrum of LSCO 
and the Fermi surface in the doping range from p = 0.03 
to p = 0.30 using the model (Al). Fixing he magnitude 
of t to 0.25 eV and the ratio t" /t' = 0.5 for all dopings 
they used p and t' as fitting parameters. Their results 
are shown in Fig. 13 The relation between p and t' as 
a function of the doping is linear to a very good approx- 
imation, 



(i(p) = -1.85* + 7.26 *'(p). 



(A2) 



Yoshida et al identified a Lifshitz transition at p « 
0.18. Our hopping parameter renormalization scheme, 
Eq. (A2), locates the Lifshitz transition at p « 0.24. 



Changing t" jt' to t" jt' = 0.25 shifts the Lifshitz transi- 
tion within our model to the experimental doping level. 
Note, that a second ARPES study by Ino et al. [24] also 
confirms the presence of a Lifschitz transition at a doping 
level of about p « 0.20. 



Appendix B: Strongly renormalized quasiparticle 
interactions 



1. Renormalization of the scattering vertex 

The calculation of the renormalized vertex in Ref. [12] 
is based on the FRG equation for the one-particle irre- 
ducible generating functional derived in Refs. [25] and 
126] which gives hierarchically coupled flow equations for 
the one-particle irreducible vertices. The RG-flow follows 
from a Wilsonian flow scheme with a sharp momentum 
cutoff. The RG equations are solved using the standard 
truncation of all vertices with more than four legs, so that 
the self energy and the four-point vertex are the only re- 
maining quantities in the calculation. For the vertex, the 
momentum dependence is discretized and constrained to 
momenta lying on the Fermi surface. The frequency de- 
pendence of the vertex is suppressed. For the RG-flow, 
the flow parameter A is taken from oo to zero. In general, 
the truncated flow tends to diverge for low temperature 
at finite energy scales A, indicating the appearance of an 
instability. The energy scale at which the divergence oc- 
curs is related to the critical temperature T c of the corre- 
sponding instability, in our case superconductivity. Note, 



however, that due to the truncation of the RG equations 
the T c is largely overestimated. 

For a study of normal state transport, the divergence 
of the vertex should be avoid. While in the experiment 
superconductivity is suppressed by a magnetic field which 
is difficult to incorporate in the RG analysis, we keep 
thge cutoff energy scale A c > finite and so suppress the 
leading d-wave pairing instability in the RG flow. The 
choice for A c is discussed in App. |B 2| 



The quasiparticle scattering vertex V r (p 1 ,p 2 ?P3?P4) = 
(P1P2 \V\P3P4 ) i s peaked for momenta p i lying close 
to the saddle-points (0, 7r) and (ir, 0) and has especially 
strong contributions for momentum transfer of (7r,7r), cf. 
Fig. 12 Further discussions of the RG flow and inter- 
pretations of vertex diagrams as in Fig. 12 are given in 
Refs. [25H2H]. 

The FRG method is a weak coupling analysis. On the 
underdoped side, below a temperature scale T* a pseudo- 
gap opens as an indication that the system is driven into 
a strong coupling phase. The strong-correlation physics 
are not included in our analysis. Hence, we restrict our 
study to dopings from optimal to strongly overdoped. 



2. Choice of parameters 



The onsite repulsion U for the cuprates is a large en- 
ergy scale of the order of the band width as at half-filling 
a Mott insulating phase is realized. Using such large val- 
ues of U, however, leads to a divergence of the RG flow 
within the normal state temperature range, since the in- 
stability temperatures are overestimated. This urges us 
to start with moderate value of the onsite repulsion of 
U ~ 2.0 1 which then is renormalized to much larger val- 
ues for the low-energy quasiparticles after integrating out 
the high-energy states. Moreover, we suppress the d-wave 
pairing instability by keeping the cutoff A c > finite in 
the RG flow, as mentioned above. For this purpose we 
choose A c = /c£T min corresponding to the lowest temper- 
ature in our study. 



Appendix C: Semiclassical treatment of normal state 
charge transport of LSCO 



1. The Boltzmann equation and its discretized 
solution 



The Boltzmann equation relates the substantial deriva- 
tive of the quasiparticle distribution function to the col- 
lision integral, taking all scattering events into account 
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FIG. 12. The renormalized vertex V(p 1 ,p 2 ,p 3 ,p 4 ) as a function of the momenta p 1 (horizontal- axis) and p 2 (vertical-axis) 
for fixed p 3 in the first angular patch. All three momenta lie on the Fermi surface and are therefore uniquely determined by 
their angular patch label. The fourth momentum variable p 4 is fixed by momentum conservation and does in general not lie 
close to the Fermi surface. We demonstrate the doping and temperature dependence by showing the vertex for three different 
values of the hole-doping p for low respectively high temperature. The initial value of the onsite repul sion is U = 2.0 1 and the 
RG-cutoff scale A c = 0.01 1. The angular discretization scheme (here 40 patches) is illustrated in Fig. 13 (a). 




FIG. 13. (a) Sketch of our curved-patch discretization of first quarter of the Brillouin zone. Plotted in color-code is the 
/o(p)(l — /o(p)), a measure of the scattering phase space. The lines represent the edges of the finite elements (boxes). For 
illustration purposes the resolution of the mesh is much lower than used for our calculations, (b) The points represent the 
relation between \i and the next-nearest neighbor hopping t' observed in the experiment [18] . We use the linear fit (red line) to 
this data as the renormalization scheme for the hopping parameters in our model, (c) Fermi surface of the tight binding model 
JaTI) with hopping parameters chosen to match the measured parameters for LSCO for different hole-doping values, cf. App. 



which yield a change of the distribution function, 
df(r,p,t) df(r,p,t) 



dt 



dt 



v(p) • V r /(r,p,i) 



-[V r V(r)]-V p f(r,p,t) 
df(r,p,t) 



dt 



(CI) 



coll 



We consider two contributions in the collision integral, 
two-particle collisions labeled by "ee" and impurity scat- 



tering labeled by " imp" : 



df(p) 



dt 



df(p) 



coll 



dt 



df( P ) 



dt 



The first contribution is expressed as, 



df( Pl ) 



dt 



(C2) 



J dp 2 dp 3 d Pi I^ 2) 3,4 *!,2,3,4, (C3) 
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with 
F 



1,2,3,4 



= \ (PlP2 \V\P3P4 ) - (PlP2 l^|P4P3 ) 
X *(Pi +P 2 -P3 -P4) 

x -jr «(e(Pi) + e(p 2 ) - e(P 3 ) " ^(pJ), ( C4 ) 



(a). For the calculations that we present here we have 
chosen a discretization of 24 patches in radial direction 
and 120 patches in angular direction. The latter are not 
distributed equally but are much finer in anti-nodal than 
in nodal direction to take care of the flat dispersion in 
the proximity of the saddle points. 



and 



1,2,3,4 



/( Pi )/(p 2 )(i-/(p 3 ))(i-/(p 4 )) 

-(l-/(Pl))(l-/(P2))/(p 3 )/(P4) 



(C5) 



where the transition rates T are generated through 
Fermi's golden rule from the matrix elements of the quasi- 
particle interaction vertex V, cf. App. [b] The second 
contribution to the collision integral is simpler due to its 
single particle form, 



df( Pl ) 



dt 



= - f dp 2 r™ p F li2 , 



(C6) 



with 

-pimp 

1 1,2 — 

and 



9tt 

\( Pl \W imp |p 2 )| 2 - 5 (s( Pl )-s(p 2 )) 



(C7) 



Fi,2 = /(Px)(l - /(p 2 )) - (1 - /(Pi))/(P2). (C8) 

We ignore electron-phonon contributions here because, 
in the cuprates, the momentum relaxation of the quasi- 
particles is dominated by electron-electron interactions. 
It is important to notice that two-particle scattering can 
only yield a finite resistivity, if momentum is transferred 
to the lattice, which only can be realized through umk- 
lapp scattering (see Sec. III). Umklapp scattering has 



special geometrical constraints and, therefore, introduces 
strong anisotropy in the Brillouin zone and requires a de- 
tailed analysis of the collision integral beyond a single- 
relaxation-time approach. 

We linearize the Boltzmann equation in terms of the 
deviation from equilibrium 5f. The collision integral is 
interpreted as a linear integral operator acting on Sf. 
We evaluate the integral kernel on discrete patches of a 
finite mesh (see Fig 13 (a)). Note in this context that 
that the energy conservation appearing as delta function 
in the scattering rates has to be treated carefully. In 
discretized momentum space, the linearized Boltzmann 
equation is reduced to a set of linear equations which are 
solved numerically. 

In order to match the Fermi surface, we use polar an- 
gle and the energy to represent each momentum vector. 
For our curved-patch discretization we can adjust the 
grid spacing in radial (energy) direction to the tempera- 
ture scale and always cover the relevant scattering phase 
space. A sketch of this discretization is given in Fig. 13 



2. Calculation of the Seebeck coefficient within 
Boltzmann transport theory 

For a study of the Seebeck effect, the Boltzmann equa- 
tion is extended by a new drift term that accounts for the 
response of the distribution function to the external ther- 
mal gradient, 



df 
de(p) 



v(p) ■ V r T 



£ = 



df(p) 



dt 



(C9) 



coll 



with the electro-chemical potential defined as S = 
V(e</> + /i). The collision integral remains unchanged. 

The distribution function that solves this equation car- 
ries a heat current and an electrical current. These two 
currents are defend as the linear response to the external 
fields as 



J i — , 



with J = 



Je 

3th 



-V r T y 
(CIO) 



Considering an open circuit geometry with j e = 0, the 
typical setup for thermopower measurement, the See- 
beck coefficient is defined as the proportionality factor 
between the applied thermal gradient and the induced 
electrochemical potential, Q = £/VT = K12/K41. For 
simplicity we assume all forces aligned along one direction 
and, thus, we can use scalar quantities. It is important to 
note that we do not make use of analytic simplifications 
as Mott's formula or Sommerfeld expansion. 



Appendix D: Conversion of computational to 
experimental units 

According to the ARPES study [18] used above to find 
the effective tight-binding parameters of LSCO, a value of 
0.25 eV for the nearest neighbor hopping parameter t fits 
the ARPES data best. This defines an energy scale which 
translates the energy scales of temperatures in units of 
t to eV. If we want to study the temperature range up 
to room temperature of 300 K the corresponding energy 
scale is ks 300 K « 0.025 eV ^ 0.1 1. In order to cover at 
least one order of magnitude in temperatures we choose 
k B T min = 0.01 1 and /c B T max > O.H. 

Our model is based on a two-dimensional interacting 
electron gas on a lattice. The experimental setup is 
based on the measurement of the resistivity of a three- 
dimensional but layered sample. For a quantitative com- 
parison of the experimental and numerical data, a con- 
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version of the two-dimensional resistivity to a three- 
dimensional resistivity is required. For a layered system 
this is very easily achieved by multiplication of the two- 
dimensional resistivity with the interlayer distance. The 



c-axis lattice constant of LSCO is given byc= 13.4 A [29] 
and the interlayer distance by half the lattice constant, 
because the unit cell of the 2-1-4-compounds contains 2 
layers. 
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